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Abstract 

A new formulation is presented to treat thermal fluctuations around the kaon 
condensate, based on chiral symmetry. Separating the zero mode from the 
beginning we perform the imaginary-time path integral to the one loop to get 
the thermodynamic potential at finite temperature. The role of the Goldstone 
mode in the kaon condensed phase is elucidated in relation to the equation of 
state. 
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Kaon condensation has been extensively studied for years JT| . As their implications there 
have been suggested many phenomena in relation to neutron stars. Among them the low- 
mass black-hole scenario, proposed by Brown and Bethe [Q], gives one of most interesting 
possibilities for explaining the current observation on the mass of neutron stars or the would- 
be neutron star in the SN1987A. According to this scenario a newly formed hot neutron star 
after the supernova explosion may collapse to a black hole (the delayed collapse) during the 
deleptonization or the initial cooling era in the star's first ~ 20s due to the softening of 
the equation of state (EOS). A numerical simulation has been already done following this 
scenario by Baumgarte et al. Through the studies of EOS it has been found that kaon 
condensation results in the large softening of EOS and a low maximum mass for neutron 
stars. It may widely lie from 1.4M to 2.1M in a relativistic calculation by Fujii et al. [|J. 
A nonrelativistic calculation by Thorsson et al also suggested similar values p[. Anyway 
kaon condensation gives a possibilty of a low maximum mass of ~ 1.4M , and thereby the 
low mass black holes can be expected around 1.4M 0. 

EOS for the kaon condensed phase has been studied mostly for low temperature. How- 
ever, the delayed collapse may occur in the hot protoneutron-star regime at temperature, 
T ~ several tens of MeV, as suggested in the simulations 0,|6[]. Hence the studies of kaon 
condensation at finite temperature is required. Since there was no EOS of the kaon con- 
densate at finite temperature, Baumgarte et al. used EOS at T ~ in their numerical 
simulation ||. Recently, Prakash et al. treated the kaon- condensed phase at finite temper- 
ature and discussed the properties of protoneutron stars within the meson exchange model, 
since there is no consistent theory based on chiral symmetry [0. 

Chiral symmetry is an important concept for kaon condensation; the kaon condensed 
phase can be represented as a chirally-transformed state from the vacuum in terms of chiral 
symmetry, and the essential results can be obtained in a model-independent way || . So we 
wish a consistent treatment of the thermal fluctuations based on chiral symmetry even at 
finite temperature. In a recent paper Thorsson and Ellis have tried to include quantum or 
thermal fluctuations within a chiral theory || . However, they have discussed only the effect 
of the zero-point energy of kaons (quantum fluctuations). Furthermore their expression for 
the thermodynamic potential looks too complicated to be tractable. 

In this Letter we present a new formalism to treat thermal fluctuations on the basis of 
chiral symmetry and demonstrate how our formalism works by calculating EOS. We shall see 
that our formalism makes the analysis of the structure about the kaon dynamics possible and 
clarifies the physics included in the condensed phase. We consider here the isothermal system 
without trapping of neutrinos for simplicity. We treat nucleons as a kind of background for 
kaons and do not take into account their dynamical degrees of freedom except the kinetic 
energy, which may follow the spirit by the heavy-baryon chiral-perturbation theory ||. 
Relativistic effects for nucleons are beyond the scope of this Letter. 

We start with the Kaplan- Nelson Lagrangian with the Goldstone field U and the Baryon 
octet B as an effective chiral Lagrangian C c hi ra i 0H3I' which consists of the symmetric part, 



C = ^ti(d^d^U) + itr{B(p - m B )B} + Dtr(B llll5 {A», B}) + Fti(B %ll5 [A», B}). 

(1) 

and the breaking part, 
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C SB = atrm q (U + U ] - 2) 

+ aitTB{£m q £ + h.c.)B + a 2 tiBB(£m q £ + /i.e.) + a 3 {tiBB}ti(m q U + h.c.) (2) 

with the quark mass-matrix, rh q ~ diag(0, 0, m s ). The coefficients a, ai measure the strength 
of the explicitly symmetry breaking: the kaon mass is given as m\ ~ 2m s a/f 2 and the KN 
sigma terms as = —m s (ai + a 2 + 203) and = —m s (a 2 + 2a 3 ). The transformation 



properties of the fields U and B under SU(3)i x SU(3)r are found in ref. |T0" |. 

It is well-known that chemical potentials, which is needed to ensure various conserva- 
tion laws in the ground state, can be introduced as the artificial "gauge fields". Here we 
must take into account the conservation laws of two quantities: electromagnetic charge and 
baryon number. Accordingly we replace the time-derivatives by the "covariant derivatives" ; 
V^qU = dU/dt + ifj,K[T e m, U],V^ =Q B = dB/dt — ifi n B + ifj,K[T e m, B], with the chemical 
potentials /iK-,//„, an d the charge operator, T em = T 3 + 1/\^3T 8 = diag(2/3, — 1/3, — 1/3) 
(T a : generators of SU(3)). Thus the partition function in the imaginary-time formalism can 
be represented as follows (r = it, (3 = 1/T); 



J chiral 



N J [dU] [dB] [dB] exp J dr J d 3 x{C chiral + 8C) 



(3) 



Here we can see that there appears the additional term in the effective action, 5C, which is 
in general non-invariant under the chiral transformation, 

6C = -£MtT{[T m ,U\fr + U[T em ,tf}} - ^tr{5t[(et[T em ,e] +e[T em ,e f ]),5]} 

x2 2 

- U^L t r{{T em , U] [T em , W]} + fi n tr{B^B} - fi K ti{B^T em , B]}, (4) 
where U = dll/dr. 

Fluctuations around the condensate are introduced as follows. First consider the chiral 
transformations on the SU(3)l X SU(3)r/ SU(3)v — 577(3) manifold. Then the eight Gold- 
stone fields (p a span the coordinates on the manifold [[□]], U(<f> a ) = exp(2iT a <p a / f) G SU(3), 
and accordingly the kaon condensed phase corresponds to the point on the manifold specified 
by Uk = exp[2z{T 4 (0 4 ) + T 5 (0 5 )}//]. [] On the other hand we know that it is attained by 
the chiral transformation from the vacuum, Uy = U(<j) a = 0) = 1, 

U K ((6}) = (U v ( = ( 2 (5) 

where £ = exp(z(M)/y / 2/) with the kaon matrix 



M 



or 


K~ 



= (0 4 ±^s)/v2 and 6 2 = 2K + R-/f 2 



Thorsson and Ellis introduced the kaon fluctuation fields in the usual manner ||, 
K ± = (K^) + K . Here we introduce the fluctuation fields by extending the idea of chiral 



x The symbol (O) implies the thermal average of an operator O. 
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rotation mentioned above. First, we generate fluctuations a around the vacuum by a chiral 
transformation Uf, Uf = rjUyr] = r] 2 with r\ = exp(iT a (p a / /). Then an additional chiral 
rotation by ( transforms Uf into the form, 

U((9)Aa)=CUf(<Pa)C- (6) 

This is our ansatz for the form of U on the manifold. It is obvious that in the limit <p a — > 
or (0) — ► 0, the matrix U is reduced to the relevant form. We can regard this procedure as 



a separation of the "zero-mode" from the full SU(3) matrix U [ 12,14] . Accordingly, the £ 
operator, which is defined by U = £ 2 , can be obtained by solving the subsidiary equation, 

£ = = UT]C, (7) 

where the matrix u is defined by the second equality in Eq.(^). 

Defining a new baryon field B' by the use of the matrix u, B' = u'Bu, we can see that 
the chiral-invariant pieces of the Lagrangian are not changed in form and only the symmetry 
breaking pieces play crucial roles, 

C(U, B) = Cq(U, B) + C SB (U, B) — > Co(U f , B') + £ 5 b(C^/C, uB'v!) 

SC(U, B) — ► 5C((U f (, uB'v)). (8) 

Thus all the dynamics among the condensate, fluctuations and baryons are completely pre- 
scribed by the non-invariant terms under the chiral transformation, Csb and SC. 

To proceed we have to find the SU(3) matrix u by Eq.(|7]). It depends on <f> a , (9) in a 
complicated nonlinear way, and thereby it is very difficult to find u for the general form of 
Uf [0. However, if we restrict ourselves to the fluctuations of kaon sector s.t. 



then we can simply find 



Uf = exp(zv / 2M//), (9) 



u — diag(/c*/|/c|, 1, k/\k\), (10) 



with 



k = oos«e>/2) cos(0/2) - |( ^jy +| sin((0)/2) sm(9/2). (11) 
For the integration measure in Eq.(0), it is invariant under the transformation (0) by £, 



[dU] — > [dUf], and further approximated by the one with flat curvature |I2]. Thus we find 



■'chiral 



f f[[dUdB'][dB'] eMS e lL(C U f , B% (12) 

i=l 



In the following calculation we only consider the kaon dynamics in nuclear matter (B — ► 
ip) in the one-loop order. As mentioned before, since we discard the dynamical degrees of 
freedom of nucleons in this Letter, we can integrate out the kaon field and the nucleon field 
separately in Eq. (O). First we consider the integration over the kaon field. Using Eq. (JlOf) 
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and retaining only the quadratic terms with respect to the kaon field in Eq. @ , we can 
write the effective action as 

Schiral = S c + S R (13) 

apart from the nucleon terms, with the tree-level contribution, 

1 



S C = (5V 



2 /4/ 2 sin 2 (£> - f(m\ -a- 2/i K b)(l - cos(6}) 



(14) 



and the fluctuation part, 

S e J f = J o drj d 3 x[- [K- - fi K cos(6)K-) [K + + fi K cos(6)K + ) - VK + VK~ 

-cos(6)(m 2 K - a)K + K- - ^ sin 2 (^)(K + + K~) 2 
- b [K + (k- - n K cos(e)K-) - K-(K + + fi K cos(6)K + )} + 0(|i^| 4 )] , (15) 

where b stands for the V-spin density, b = (V>7og(3 + r 3 )ip) / f 2 and a the KN sigma-term 

contribution, a = (M)^ KN ijj) / ' f 2 with E KN = I „° ). R Note that Eqs. (14), (15) 

are the exact consequences from the general argument given in Eq. (8). Expanding the 
fluctuation field with the Matsubara frequency u n = 2imT and performing the integration, 
we find the kaon-loop contribution, 

\nZ e K ff = Y,H\D eff \- l/2 ), (16) 

n.p 

with the inverse thermal Green function D e ^ , 

D eff = 2 ( + ^k){^+ + ^k)+^k sin 2 (6) 2(jl K + b)u n \ . , 

P \ -2{fi K + b)u n u 2 n + (cD_ -p. K ){u>+ + p, K ) ) ' 1 ; 

where Q± = ±6+ (p 2 +fh*^ +b 2 ) 1 / 2 , Jik = fi>K cos{6) and the effective mass, m*^ = cos(9)[m 2 K — 
° ■ □ 

The excitation energy of the kaonic modes is given by the poles of the thermal Green 
function. Then we have two solutions corresponding to the modes, 



El = (ci + c 3 + 2c 2 ) ± J{c 3 + 2c 2 ) 2 + 4 Cl c 2 , (18) 



2 It would be interesting to compare our result with the one given in ref . § ; although our expression 
might look rather simple, the results, e.g. the dispersion relations for the kaonic modes, are 
consistent with them. 

3 In the limit (9) — ► Q,u>± — > oj± = ±6 + (p 2 + m*^ + b 2 ) 1 / 2 , which is nothing but the dispersion 
relation for kaons in the normal phase [|],||] . 
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where c\ = (&_ — C2 = j^K + b and C3 = l/2/i|- sin 2 (6 ) ). In the condensed phase 
((#) 7^ 0), the mode corresponding to E_ is the Goldstone mode as a consequence of the 
breakdown of V-spin symmerty in the ground state; the ground state is not an eigenstate of 
the conserved operator V3. It is easy to show E_(p = 0) ~ 0. In the tree-level approximation, 
c\ = as a result of the extremum condition for S c , dS c /d(9) = 0, which means £L = 0, 
while we should find E_(p = 0) = + 0(h) once taking into account the fluctuation effects 
(kaon loops) ||. When we consider the thermal kaon loops, E_ directly enters into the 
Bose-Einstein distribution function and it should diverge at p = 0. The other is the massive 
mode and can be neglected for temperature we are interested in (T < lOOMeV). 
Moreover, if we can safely neglect the C3 term, we get a simple expression for E±, 



E±=ti±±fi K = ^/p 2 + m* 2 + b 2 ±(b + jl K ). (19) 

We can expect the pertinence of this approximated formula qualitatively by several reasons: 
first, this term includes the KK scattering amplitude in the leading order in difference 
from the other terms in Eq. (|17|), which should be small at low-energy regime due to the 
Goldstone-particle nature of kaons. Secondly, in the limit, (9) — > or fi K — ► 0, c 3 gives 
no contribution. We also know that (9) and fix are inversely proportional to each other 
referring the previous studies, e.g., Leaving the detailed discussions in ref. |TJ|], we 

find that the approximation is justified if the following condition is satisfied, 

c 3 « ^TC 2 (3/2) or ^sin(^)«i(27rCT) 1 / 2 C(3/2), (20) 

with C 2 = m*x + b 2 . We shall see, by numerical calculations, that this approximation really 
works well even at finite temperature except very low temperature, where thermal effects 
should be trivially unimportant. Finally, it is to be noted that this approximation never 
violate the Goldstone-particle nature, £L(p = 0) = + 0(h). 

Then the thermodynamic potential by the kaon-loop contribution, Vt e ^ = — TlriZ^, 
reads 



V I Wr^ [E ^ ) + E -^ ]+TV I (0alB(l- C -^-W)(l- C -^-W). (21) 



The first term is the zero-point energy (ZP) contribution and the second one the thermal 
one. Since ZP contribution includes the divergent integrals, they should be renormalized in 
some way. Following the method proposed in ref. ||, we can extract the finite contribution 
Q zp,r ; under the approximation, C3 = 0, Q zp,r takes a simple form, 



n zp < r l — v 

64tt 2 



(m 2 K - 3C 2 )(m 2 K - C 2 ) + 2C 4 ln 



C 2 



(22) 



We shall see numerically that Q zp,r results in a tiny value, as consistent with ref. ||, we 
discard it in the following. 

In the rest of this Letter we present some results about the dispersion relation and EOS at 
finite temperature. For nucleons we should take into account at least two effects, otherwise 
the results should become meaningless. First one is the thermal effects for the kinetic energy 
term; it is given by the standard formula, 
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n nn „ _ 2TV ^ J + e-K*-K)), (23) 



where e = p 2 /2M for nonrelativistic nucleons and p p is defined by the relation, p p = p n —p K . 
Secondly, it is well-known that the nuclear symmetry energy plays an important role for the 
ground-state properties of the condensed phase. Since we have already included the kinetic 
energy for nucleons, we take into account only the potential energy contribution. Following 



Prakash et al. [13| we effectively introduce a symmetry energy contribution, 

n sy mm = VpB (i_ 2x) 2 S pot (u), (24) 

where ps the baryon number density, u = Pb/Po (Po — 0.16fm~ 3 the nuclear density), and 
x is the proton-number fraction. The coefficient S pot (u) is given as S pot (u) = (So — (2 2 / 3 — 
l)(3/5)Ep)F(u) with the empirical symmetry energy Sq ~ 30MeV and the Fermi energy 
Ep at pQ. We take here F(u) = u as a form of the control function for simplicity. Thus the 
nucleon contribution is given by Vt e ^ = Vt k ™ + Q s ^ mm . 

Finally, the total thermodynamic potential Qtotai is given by further adding the one for 
leptons (electrons and muons), Qi, 

^total — + ^//^ + ^tv^ + (25) 



where Q c is the tree-level contribution, Q c = —TS C , given by Eq. ( |14"D and we can use 



the formulae of free leptons for Qi. The parameters x, {6) and the chemical potential px 
are determined by the extremum conditions for given density and temperature like in the 
T = case, e.g., []l]|||5||: dVt to tai/dx = 0, dtt to tai/d(6) = and dVt to tai/dp K = by the use of 
Eq.©. 

In the following we present some numerical examples to see how our formalism works 
and how large the thermal effects are, leaving the full discussion about the thermal effects 



and the parameter dependence in a separated paper fllj]. In numerical calculations we use 
the values, aim, = -67 MeV, a 2 m s = 134 MeV and a 3 m s = -222 MeV §. First of all 
we examine the pertinence of the approximation C3 = 0. We compared the exact dispersion 
relation (|i~8D with the approximated one (|19| ) as a function of momentum up to 500 MeV 
and checked that this is a fairly good approximation. In Fig.l we show a typical example, 
which might be the worst case at T = 50 MeV in such a sense that the difference between 
them becomes largest. We can see that the difference is very small; in fact, it causes the error 
of only less than several MeV-fm -3 for pressure. As a general trend the difference becomes 
smaller as density or temperature becomes higher (c.f. Eq. ( p0"D ). It should be worth 
noting that once this approximation works, the expressions of thermodynamic quantities 
like charge, entropy, energy become very clear [O . 



In Fig.2 we present an example of EOS for the hot neutron-star matter by using C3 = 0. 
We construct EOS for the condensed phase, referring to the one for normal neutron-star 
matter; P = Py + P c hirai- The pressure P c hirai is simply given by Eq. (|25|) through the 



4 We also compared our dispersion relations with ref. ||] at T = and found that they are consistent 
with each other regardless of the magnitude of (6) in the ground state. 
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thermodynamic relation, P c hira\V = —Qtot- As a reference EOS for normal neutron-star 
matter we employ here the one suggested by Prakash et al. (T^|; Py = u 2 d(V(u)/u)/du with 



the potential contribution V(u). The results indicate the first-order phase transition in Fig. 
2. Hence, the Maxwell construction should be imposed to get the real EOS. We can easily 
see that the thermal effects are remarkably larger than those in the normal phase; this is 
caused by the reduction of the chiral angle, and the soft mode mainly contributes to this 
effect. 

Finally we have checked the smallness of the ZP contribution to be tiny; the value of the 
ZP contribution to pressure, Pzp = —Q zp ' r /V, is always less than several MeV-fm -3 over 
the density-temperature region we are interested in. 

In this Letter we have presented a formalism to treat fluctuations around the con- 
densate within the framework of chiral symmetry. Our approach is based on the group 
theoretical argument; Goldstone fields are regarded as the coordinates on the chiral 
SU(3)i x SU(3)r/SU(3)v manifold and the kaon condensed phase can be represented as 
a chiral rotated state from the vacuum. Using this idea we have introduced fluctuations 
around the condensate by way of the successive chiral transformations. We derived an effec- 
tive action and extracted the dispersion relation for kaonic excitations. There appear two 
modes: Goldstone-like soft mode and very massive one, which correspond to K~ and K + 
mesonic excitations, respectively. We have discussed that the form of the dispersion rela- 
tion can be reduced to a simple one by the approximation, C3 = 0, which was also checked 
numerically and proved to be fairly good. 

Applying our formalism to the derivation of EOS as an example, we have pointed out 
the largeness of the thermal effects in the condensed phase; they result in a remarkable 
reduction of pressure in the condensed phase, compared with the normal phase. 

In this Letter we have not taken into account the dynamical degrees of freedom of 
nucleons consistently, but it is not difficult to treat them consistently based on the idea 
given in this Letter |JL4|j. 



We thank T. Muto for discussions and careful reading of this manuscript. This work was 
supported in part by the Japanese Grant-in-Aid for Scientific Research Fund of the Ministry 
of Education, Science and Culture (08640369). 
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Figure Captions 

Fig. 1 Dispersion relation for kaonic modes in the condensed phase at u — 3.87 and T = 50 
MeV. Solid lines show the results given by Eq. fliTf), while dotted lines by Eq. (|T%D. 

Fig. 2 Equation of state for neutron-star matter at finite temperature. EOS at T=0 , which 
is the same one given in ref. ||, is also shown for comparison. Dotted lines show EOS 
for normal neutron-star matter. 
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